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The ALADIN experiment aims at observing how the critical magnetic field of a superconducting 
Aluminum film is modified, when it constitutes one of the reflecting surfaces of a Casimir cavity. If 
successful, such an observation would reveal the influence of vacuum energy on the superconducting 
phase transition. In this paper a rigorous analysis of experimental data is reported, the results are 
discussed and compared with theoretical predictions based on Lifshitz theory of dispersion forces, 
and the BCS formula for the optical conductivity of superconductors. The main novelty with respect 
to a previous data analysis by some of the authors, is the use of a cross-correlation method which 
is more rigorous and leads to better estimates. 
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I. INTRODUCTION 

The Casimir effect U provides a striking manifesta- 
tion of vacuum quantum fluctuations of the electromag- 
netic field in bounded geometries, and it represents a rare 
example of a purely quantum phenomenon that can be 
tested at the mesoscopic scale. The Casimir effect has 
received much attention in the past decades, thanks to a 
wave of new experiments which made it possible to mea- 
sure the Casimir force with unprecedented precision. For 
a recent review of these experiments and a critical survey 
of the numerous theoretical investigations on the Casimir 
effect in materials, we address the reader to the recent 
monograph [3j. 

Despite the impressive theoretical and experimental 
advances made in the past twenty years, the Casimir ef- 
fect still faces important unsolved questions at the fun- 
damental level, in particular the problem of reconciling 
the vacuum energy density and its interaction with the 
gravitational field, known as the cosmological constant 
problem 0, H|. No experimental verification that vac- 
uum fluctuations gravitate according to the equivalence 
principle has been obtained so far, even though there are 
theoretical expectations that this should be the case 
[13 ]. Relying upon these considerations, some of us stud- 
ied the effect of a gravitational field on a rigid Casimir 
cavity, by computing the net force acting on it: interest- 
ingly, it was found that a Casimir apparatus, when sub- 
ject to the weak gravitational field of the Earth, should 
experience a tiny push in the upwards direction [15-18]. 
In [l5| it was argued that an experimental verification 
of this effect is extremely hard, if not impossible, un- 
der static conditions, by virtue of the extreme smallness 
of the expected force. A better possibility would be to 



carry out the measurement dynamically, i.e. by modu- 
lating the Casimir energy stored in a rigid cavity in a 
known way. Such a modulation of the Casimir energy 
could be achieved by altering periodically the reflectiv- 
ity of the plates. Recently, a significant modulation of 
the Casimir force between a highly doped semiconduct- 
ing membrane and a gold plate has been demonstrated 
experimentally [l9|, by shining periodic laser pulses on 
the membrane which determine a large change in the 
charge carrier density of the membrane. Despite being 
interesting for Casimir studies, this result is not suitable 
for "weighting" aims: the energy supplied to the system 
to induce the change in carrier density is many orders 
of magnitude larger than the variation of Casimir en- 
ergy. This would make it extremely difficult to observe 
the tiny fraction of mass change due to the Casimir con- 
tribution. On the contrary, in our scheme based on the 
superconducting phase transition, the total change of en- 
ergy is of the same order of magnitude as the change in 
the Casimir energy and therefore its contribution might 
in principle be observed. 

This is the framework of the ALADIN experiment, 
whose aim is to observe the variation of the Casimir en- 
ergy stored in a superconducting Casimir cavity, consti- 
tuted by a thin superconducting film separated by a thin 
oxide layer from a thick gold substrate, across the super- 
conducting phase transition. The scheme of detection is 
based on a measurement of the critical magnetic field that 
destroys the superconductivity of the film, whose magni- 
tude is expected to be influenced by the Casimir energy. 
If successful, the experiment would thus reveal the influ- 
ence of vacuum energy on a phase transition. Another 
distinctive feature of our setup is that we use rigid cav- 
ities, that are obtained by deposition techniques, a fea- 
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ture which might be useful to investigate experimentally 
the dependence of the Casimir energy on the geometrical 
shape of the intervening bodies, an issue that remains 
under scientific debate also at a theoretical level Q. 

The plan of the paper is a follows: in Sec. II we briefly 
describe the experimental setup and the measurement 
method, Sec. Ill is devoted to the analysis method based 
on cross-correlation. Finally, the experimental results are 
discussed in Sec. IV. Indeed, a preliminary analysis of the 
data reported in this paper has already been performed 
some time ago, but without the necessary rigor, and we 
reported on it in [20]. The analysis reported in this pa- 
per, based on a procedure of cross-correlation, is more 
rigorous and the results are now better estimated. 



II. ALADIN: EXPERIMENTAL SETUP AND 
EXPECTED EFFECT 

Before we describe our experimental setup, it is useful 
to briefly recall the principle at the basis of our experi- 
ment, that was described in detail in the works [2l|, l22j . 
The starting observation is that according to Lifshitz the- 
ory [3|, the Casimir energy is determined by the opti- 
cal properties of the plates. Since the optical proper- 
ties of a superconductor are sharply different from those 
of a normal metal [23|, one is led to expect that the 
Casimir free energy stored in a superconducting 

cavity should change across the superconducting tran- 
sition. The change AF^ = F^ C) - F S (C) of Casimir 
energy was estimated in [2l|, [22| (see also 0, [HI), on 
the basis of Lifshitz theory by using the BCS formula 
for the optical conductivity of superconductors, and it 
was found to be extremely small. This is not surprising 
of course, because the superconducting transition alters 
the optical properties of a metal only in the microwave 
region, which constitutes a small window in the wide fre- 
quency range that contributes to the Casimir energy. The 
latter typically extends up-to a few times characteristic 
cavity frequency uo c = c/2d, with d the plate separa- 
tion, which for typical submicron separations belongs to 
the infrared region of the spectrum. The smallness of 
the fractional change of Casimir energy across the super- 
conducting transition makes it impossible to observe the 
corresponding change in the Casimir force on the plates, 
with present day sensitivities in force measurements. 

The Aladin experiment uses a detection scheme which 
does not involve at all a force measurement, as it aims 
at observing how the variation AF^ of Casimir energy 
influences the critical magnetic field H c of a thin super- 
conducting film which is part of a Casimir cavity. To see 
how this comes about, we recall [26| that the magnitude 
of the parallel critical field H c for a thick superconducting 
slab of volume V can be determined by equating the mag- 
netic work VH^ I 8tt required to expel the magnetic field 
from the sample, with the so-called condensation energy 
^cond(^) of the material, which represents the difference 
of Helmholtz free energies between the normal and the 



superconducting phases: 
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(2.1) 



When the film is one of the two plates of a Casimir 
cavity, we have to augment the right-hand side of the 
above Equation by the difference AF^ between the 
Casimir energies in the normal and in the superconduct- 
ing phases: 
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(2.2) 



In writing this relation, we are tacitly assuming that 
the fluctuating electromagnetic field in the Casimir cav- 
ity does not alter significantly the properties of the su- 
perconductor, and in particular its condensation energy. 
This may be considered as a plausible assumption, as 
far as AF^ C \T) <C e cond (T). According to Eq. (|2^|) . 
the variation AF^ C \T) of Casimir energy determines a 
change in the magnitude of the critical field H C1 which 
for AF^ C \T) <C £ CO nd(T) is estimated to be: 



SH C _ AF(°\T) 
H c 2£ cond (T) 



(2.3) 



The key thing to notice is that the condensation energy of 
a thin superconducting film can easily be orders of mag- 

(C) 

nitudes smaller than typical Casimir energies Fn , and 
therefore one may hope that even tiny fractional changes 
of Casimir energy AF^ C \T) can determine observable 
shifts of the critical field. For example, for a Beryllium 
film, we estimated [U, [22| that a relative variation of 
F(°) of one part over 10 8 might lead to a 5 percent vari- 
ation of critical magnetic field. 

Detailed numerical computations [22| show that the 
magnitude of the effect increases for thin films, because 
they have a smaller condensation energy, and for small 
cavity widths d, because the change of Casimir energy 
becomes larger. It is important to remark that the rel- 
ative shift of critical field was found to be roughly pro- 
portional to the inverse of the transition temperature T c 
of the superconducting material. The main reason why 
superconductors with low T c lead to a larger effect is that 
the condensation energy is empirically known to scale as 
T c 2 6 [27], and this makes the denominator in the r.h.s. of 
Eq. (|2.3p decrease faster than the numerator as T c de- 
creases. As a compromise between the possibility to per- 
form preliminary tests, easy change of structures, statis- 
tics and signal-to- noise ratio, we chose to work around 1.5 
K, using Al as superconducting material. Besides having 
a low critical temperature, Al is a material that oxides 
easily, and this makes it easier to grow oxide layers of 
controllable thicknesses, well attached to the supercon- 
ducting film, that constitute the dielectric medium of our 
metallic Casimir cavities. The configuration used in the 
experiment is a three-layer cavity, made of a thin super- 
conducting Al film (5 ~t~ 10 nm), a thin dielectric layer 
of native oxide (A^Os)^ 10 nm), and a thick metallic 
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layer of Au (100 nm). However, at present we are consid- 
ering different configurations in order to obtain a larger 
signal-to-noise ratio in the expected effect. 

In our experiment, we used a cryogenic system based 
on the Heliox VL 3 He cryostat, inserted into a dewar 
equipped with magnetic screening, which isolates the 
samples from external EM fields. Since it is extremely 
difficult to keep the cryostat temperature perfectly 
constant, we did not try to measure the critical field 
H c of the samples as a function of the temperature. 
Rather, we measure how the critical temperature 
T C (H) of the samples changes as a function of the 
applied magnetic field H. More precisely, for each value 
of the applied field H we measure the relative shift 
St = (T C (H) — T c o)/T c o in the critical temperature T c 
of the sample, with respect to the critical temperature 
T c o of the same sample in zero field. As we shall explain 
below, by using the method of cross- correlations the 
shift St can be measured much more accurately than 
the individual critical temperatures in the applied and 
in the null fields. Our theory predicts (see Figure [T] 
and comments below) that the curve H(St) for the bare 
film should lie below that for the film in the Casimir 
cavity. Since this is a differential measurement, we need 
a very good sensitivity in temperature, of order a few 
fiK in the case of Al. As described in [20j |. the critical 
temperature is determined by measuring the resistance 
of a sample R(T) in a four-wire configuration around 
the phase transition, for different external applied fields. 
Several measuremements have been performed and 
different samples have been tested, so as to find the best 
experimental conditions for a good signal-to- noise ratio. 
For a good data analysis to be possible, it was important 
to make sure that the transition curves did not change 
their profiles in time, or after applying a magnetic field. 
By employing these criteria, the best samples have been 
selected. Data reported in the following were obtained 
from the samples showing the sharper transition and 
the highest homogeneity among transition curves. For 
a detailed description of sample preparation, cryogenic 
apparatus and measurement scheme, see [20|. 

The expected effect is shown in Fig. [TJ Here the critical 
magnetic field H is plotted against the shift St = 1 — t of 
the reduced critical temperature t = T C (H) /T c q (with T c o 
zero-field transition temperature). We recall [26j that for 
a (bare) thin superconducting film, with a thickness D 
much less than the penetration length A, incomplete field 
expulsion leads to higher values of the parallel critical 
field H, as compared to bulk samples. For 1 — t <C 1 
(as in our case), the thin- film critical magnetic field H 
follows the simple law [26| 

H = HoVM^-VT^t , (2.4) 

where A(0) is the penetration depth and Hq is the bulk 
critical field, both determined for zero temperature. 
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FIG. 1: Simulation of the expected signal for a bare thin Al 
film of thickness D = 14 nm (lower curve) and for a cavity 
consisting of a similar Al film, covered by a 6 nm dielectric 
layer and a 100 nm Au mirror (upper curve). 



For the in-cavity film we can divide the expected sig- 
nal curve (upper curve in Fig. [1]) into three temperature 
regions. For temperatures far from the transition tem- 
perature of the film (region not shown in Fig. [T]) the 
in-cavity curve coincides with the bare film curve, since 
the vacuum energy contribution becomes negligible com- 
pared to the condensation energy. When St ~ 3 x 10 -5 , 
the change of Casimir energy is small but no longer neg- 
ligible, and a perturbative approach is possible in SH/H: 
for H w 5 -=- 6 mT the two curves are expected to differ 
by an amount St ~ 6x 10 -6 . Note that, since the critical 
temperature of Al is 1.5 K, this corresponds to a shift in 
temperature of order 10 fiK. Finally, for lower tempera- 
tures the dependence of St on H for the in-cavity curve 
is expected to differ from the bare film case, because the 
Casimir energy contribution is of the same order of mag- 
nitude as the condensation energy. In this region we are 
not able to perform any perturbative calculation, and 
there is no theoretical prediction for the curve's trend 
(dashed line in Fig. [1]). 



III. ANALYSIS METHOD BASED ON 
CROSS-CORRELATION 

As was explained in the previous Section, for each value 
of the applied magnetic field H we need to determine ac- 
curately the fractional shift St of the critical temperature 
of the sample, relative to its critical temperature T c o in 
zero field. As described in Ref. [20], the major source 
of noise in our measurements is the electronic noise at 
the read-out amplifier. This noise has a "fast" compo- 
nent (with time scale of one second) that determines a 
statistical error of a few fiK, and a slow thermal drift 
(linear in time, about 50 \iK per hour) that produces a 
shift proportional to the time elapsed between two mea- 
surements. To correct for the thermal drift we arranged 




---l-> R(T, H =0;ti) 

2-> R(T,H;^+k) 

--~3-> R(T, H =0;t 2 ) 
R (T, H =0;^) 
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FIG. 2: Example of the measurement sequence. All measure- 
ments are performed at equal time intervals. Before and after 
each measurement with the applied field, a measurement in 
the absence of field is performed. From the two measurements 
in absence of field, it is possible to reconstruct the zero- field 
curve that would be measured simultaneously to the tran- 
sition curve in the presence of the magnetic field. AT C is 
the distance in temperature between the R(T) curve in the 
presence of the magnetic field (continuous line), and the re- 
constructed R(T) curve in zero field (dashed-dotted line). 



our measurements in series of triplets, as is shown in Fig. 
[21 Each triplet consists of three measurements of the 
curve R(T) that are equally spaced in time, the first and 
the last one of which (left-most and right-most dashed 
curves in Fig. [2]) are performed in zero-field, while the 
intermediate one (continuous line in Fig. [2]) is done with 
the field applied. Since the slow thermal drift is linear 
in time, and since the two zero-field measurements are 
performed at equal time intervals before and after the 
applied field measurement, it is possible to "reconstruct" 
out of them the position that the zero-field curve R(T) 
would have occupied in the R — T plane (dotted-dashed 
line in Fig. [2]), had it been measured at the same time as 
the applied field curve. The relative shift St in the crit- 
ical temperatures is then determined by comparing the 
measured curve R(T) in the applied field, with the recon- 
structed curve in zero field. Having explained the general 
scheme of the measurements, let us see how the method 
of cross-correlations permits to accurately determine St. 

The key feature of the R(T) curves that allows to use 
the cross-correlation method to accurately determine the 
shifts St is that, for each sample, their shape is prac- 
tically independent of both the intensity of the applied 
magnetic field, and the time at which the measurements 
are performed. In other words, the curves R(T) that are 
taken at different times, or in different fields, appear to 
differ from each other just by some horizontal transla- 
tion A along the temperature axis. The cross-correlation 



method is ideally well-suited to determine the amount 
of this translation, independently of any model for the 
shape of the curves. The idea is very simple, and consists 
in looking for the translation that maximizes the overlap 
between any two R(T) curves. Let us see how this works 
out in detail. In reality, each R(T) curve consists of a 
large number of data points more or less scattered in the 
T — R plane, around some ideal transition curve. We 
consider only data points that belong to the transition 
region, having a width of a few mK, around the critical 
temperature. Typically, this region contains a few thou- 
sand points for each curve. At this point, we cover the 
T — R plane by a rectangular grid whose axis are parallel 
to the T and R coordinate axis, and whose steps are St 
and sr respectively. The points pij of the grid in the 
T — R plane thus have coordinates {i St + «t, j sr + a#}, 
where {(1t,clr} are the coordinates of the grid's origin. 
We let Hi j be the number of data points that occupy 
the grid cell whose left-down corner coincides with the 
point p(i,j) of the grid. Clearly, the number of points 
occupying the non-empty cells depends on the size of the 
cells, i.e. on the steps St and sr. We shall discuss below 
the criterion we used to choose these steps. In practice 
for the chosen size of the steps, the occupied cells turn 
out to include about five or six data points. Consider 
now any two transition curves R^ (T) and R^ (T) (not 



and T-L^} be the re- 



necessarily distinct), and let n\ 3 

spective histograms. For any translation of R^ b \T) by 
n steps along the T-axis, we define the cross-correlation 
C( a ' 6 )[n], to be the quantity: 



(3.1) 



The quantity C (a) [n] = C (a ' a) [n] shall be denoted in what 
follows as the self-correlation of the curve R^ a \T). Intu- 
itively, C (a ' 6) [n] measures how well R^(T) and R^(T) 
overlap, after we translate R^ b \T) by an amount n x st 
along the T-axis. 

In Figs. [3] and 3] we plot the normalized self-correlation 
C[n] of one of our R(T) curves, for two different steps st- 
In Fig. [5] we plot the self-correlations and the crossed 
correlations of two of our curves, one measured in zero 
field and the other in a non-zero field. We verified that 
all correlations have a Gaussian shape around the central 
peak. More importantly, we see from Fig. [5] that the 
maximum of the cross-correlation between the zero-field 
curve and the non-zero-field curve is of the same order 
of magnitude as the maximum of the self-correlations of 
the two curves. From this circumstance we infer that 
the two curves actually have equal shapes and that they 
can be made to overlap by a translation AT along the 
temperature axis. From the cross-correlation plot, one 
can estimate AT = n max x s T , where n max is the shift 
for which the cross-correlation reaches its maximum. 
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FIG. 3: Self-correlation of the R(T) curve obtained with 
lattice step on the abscissa equal to 15/iK. 
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FIG. 4: Self-correlation of the R(T) curve with lattice step 
equal to 5/iK. 



A. Grid sizing and error estimates 

To choose the best size of the grid cells, we investi- 
gated the behavior of the self-correlation on a trial curve 
for several different choices of the steps St along the tem- 
perature axis. We pointed out earlier that the distribu- 
tion of bins around the autocorrelation peak is Gaussian 
shaped. We found that the width a of the Gaussian de- 
creases with the grid step until it stabilizes, as shown in 




Cross-correlation 1-2 
Self-correlation 1 
Self-correlation 2 



FIG. 6: Width a of the autocorrelation versus the grid step 



Fig. [6l On the basis of this behavior, a step of 15/liK in 
the plateau region was chosen, which is still large enough 
to ensure that the typical number of points in the occu- 
pied cells of the grid is five or six. 

The error on the temperature shifts AT between any 
two curves R(T) was estimated by using the familiar 
jackknife method. We divided the transition region of 
the T — R plane into n stripes parallel to the T-axis, 
and of equal widths along the i?-axis, and we covered 
each of these stripes with grids of equal steps st and 
sr. The data points falling in each of the n stripes 
were then analyzed by the cross-correlation method, pro- 



viding n independent estimates AT^ n \ k = l,...,n 
the temperature shifts. For each number of subdivi- 
sions n, we then determined the corresponding average 



of 



shift AT^ 



(n) 



jk=l 



and the variance a 
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AT 



En 



(Ar[ n) -ATW) 2 



To find the optimal number of 

stripes, the above process was repeated for n<G {1,40}. 

In Fig. Owe plot the behavior of versus n. As we 

See, (J y\ rji reaches a plateau for n w 10 — 15, and on this 
basis we adopted n = 15 for our final assessment of the 
errors. 

Having determined the optimal grid size and number 
of stripes, we could then determine the best estimates for 
the differences AT C = T C (H) — T c q between the critical 
temperatures in the presence of the magnetic field and 
in null field, with the relative errors. For each triplet of 
measurements as described in Sec. Ill (see Fig. [2]), we 
estimated 



ST(nK) 



AT C = AT (1 



(1,2) 



AT d,3) 



(3.2) 



FIG. 5: Self-correlation of two R(T) curves, green and blue, 
and cross-correlation among the two. 



where the superscripts 1 and 3 refer to the curves in zero 
field, and the superscript 2 to the curve in the presence 
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FIG. 7: The error cr^^ versus the number n of subdivisions. 
The error reaches a plateau for n w 10-15. 

of the applied field. The error on AT C was taken to be 



J AT C = y ^AT (1>2) + ^2^(1.3) J ~ COV (AT (1 , 2) ,AT (1 ,3)) 

(3.3) 

where (Jat 1j2 an d ^ATi^ are the standard deviations ob- 
tained as explained previously with n = 15. 

IV. EXPERIMENTAL RESULTS 

The first thing that we checked is that the data for the 
bare film actually follow the theoretical law Eq. (j2.4p . 
According to that relation, the shifts ST should have a 
parabolic dependence on the applied magnetic field H. 
This expectation is fully verified by our data, as can be 
seen from Fig. where the data are plotted together 
with a parabolic fit (continuous line). 

The data for the in-cavity film are shown in Fig. [9] 
(note the different scales for H and ST in comparison 
with Fig. [8]) together with a parabolic fit on higher mag- 
netic field data. It is apparent that low-field in-cavity 
data show deviations from the parabolic behavior, unlike 
the bare-film data. 

The different behavior of in-cavity data compared to 
bare-film data can be better appreciated from Fig. [TUJ 
where they are both plotted. It should be noted that in 
Fig. [10] the shifts ST/T C are reported as a function of 
the absolute value of the magnetic field, the shifts being 
independent of the sign of H. The two curves in Fig. 
[TUl are the same as in Fig. [TJ We observe again that 
the bare-film data lie nicely on the expected theoretical 
parabolic curve. More detailed comments are in order for 
the in-cavity data in Fig. [TDJ Let us consider first the re- 
gion corresponding to H « 5-1-6 mT: as discussed in Sec. 
II, for these larger fields the Casimir energy variation is 
sufficiently small compared to the condensation energy to 
justify our perturbative calculations. In this region one 
expects a small deviation of the in-cavity data from the 
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FIG. 8: Bare film data: the points follow the expected 
parabolic behavior. 



bare-film parabolic behavior. Unfortunately, the error 
bars are of the same order of magnitude of the expected 
small deviations, and therefore a better sensitivity would 
be needed to ascertain the effect in this region. For lower 
magnetic fields, the condensation energy and the vari- 
ation of Casimir energy become of the same order, the 
perturbation expansion is no longer possible, and devi- 
ations are possibly larger. Interestingly, in this energy 
region, the data of the in-cavity film are no longer com- 
patible with a parabolic behavior, as is the case for the 
bare film. This suggests the hypothesis that in-cavity 
data display an anomalous trend with respect to bare 
film data. Note that the preliminary analysis reported in 
[20| . even if less rigorous, is compatible with the present 
results. However, it is not possible to draw final conclu- 
sions from this observation, since error bars are too large 
and a theoretical model for low condensation energy is 
missing. 

To sum up, in-cavity data show a behavior which is 
qualitatively different from the bare film case, and com- 
patible with the predictions of Refs. [2l|, [22|. Even if we 
cannot draw definite conclusions from the present anal- 
ysis, its results encourage us in continuing this research, 
and exploring different configurations to enhance the ef- 
fect. 
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FIG. 10: Theoretical prediction and experimental results. In- 
cavity film data (squares), bare film data (diamonds). The 
lower curve shows the theoretical prediction for bare film data, 
the upper one that for in-cavity film data. The point-dashed 
line indicates the region where a definite theoretical prediction 
is not possible. 



FIG. 9: In-cavity data together with a parabolic fit on higher 
magnetic field data: for low magnetic field, corresponding 
to the dashed part of the curve, the data do not follow the 
parabolic behavior. 
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